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We develop a dynamic theory of output coupling, for fermionic atoms initially confined in a 
magnetic trap. We consider an exactly soluble one-dimensional model, with a spatially localized 
delta-type coupling between the atoms in the trap and a continuum of free-particle external modes. 
The transient dynamics of the atoms, as they leave the trap, is investigated in detail. Two impor- 
tant special cases are considered for the confinement potential: the infinite box and the harmonic 
oscillator. We establish that in both cases a bound state of the coupled system appears for any 
value of the coupling constant, implying that the trap population does not vanish in the infinite- 
time limit. For weak coupling, the infinite-time spectral distribution of the outgoing atoms exhibits 
peaks corresponding to the initially occupied energy levels in the trap; the heights of these peaks 
increase with the energy. As the coupling gets stronger, the infinite-time spectral distribution is dis- 
placed towards dressed energies of the fermions in the trap. The corresponding dressed states result 
from the coupling between the unperturbed fermionic states in the trap, mediated by the coupling 
between these states and the continuum. In the strong-coupling limit, there is a reinforcement of 
the lowest-energy dressed mode, which contributes to the spectral distribution of the outgoing beam 
more strongly than the other modes. This effect is especially pronounced for the one-dimensional 
box, which indicates that the efficiency of the mode-reinforcement mechanism depends on the steep- 
ness of the confinement potential. In this case, a quasi-monochromatic anti-bunched atomic beam 
is obtained. Results for a bosonic sample are also shown for comparison. 



PACS numbers: 03.75.Pp, 03.75.Ss, 32.80.Pj 

I. INTRODUCTION 

The demonstration of the first atom lasers [1|, 
has led to questions that are reminiscent of those asked 
when the first optical lasers were put to work. What is 
the dynamical behavior and the statistics of the outgo- 
ing beam? How monochromatic it is? Atoms offer an 
interesting twist to these questions, since they may have 
bosonic or fermionic behavior, while for photons only the 
bosonic character manifests itself. One may then ask how 
the statistical properties of the trapped atoms affects the 
outgoing beam. 

At zero temperature, one may guess that the behavior 
of the outgoing beam should be markedly different in the 
two cases, since for fermions there would be a multitude 
of populated trapping levels, due to the Pauli exclusion 
principle, while for a bosonic gas all the atoms would be 
in the ground state. One expects therefore that fermionic 
systems should exhibit a richer dynamics, at zero tem- 
perature, as compared to bosonic systems, which have 
been described by one-level models [a ItLH fgllTollTllIT^ 
or mean field theories Also, coherence proper- 

ties of fermionic beams are expected to be quite differ- 
ent from their bosonic-beam counterparts. Indeed, while 
bosonic beams coming from thermal sources exhibit a 
bunching effect, anti-bunching has already been experi- 
mentally demonstrated for electron beams [T^ ■ 

A simple model for a beam of fermionic atoms ex- 
tracted from a trap was analyzed some years ago [Tsj . 
More recently, the so called input-output formalism de- 
veloped for photons 0,1131 and applied to bosonic atoms 



[g was generalized to fermionic species |2]| . 

One should note that effects concerning the multi- 
level structure of the trap should also appear in a non- 
mean field theory for bosonic atom lasers with a non-zero 
temperature, problem that has been addressed very lit- 
tle, and only within the Popov approximation p^j for 
trapped systems jl^, 0] . 

Theoretical work on degenerate fermionic gases has 
been greatly stimulated by the first propositions of a su- 
perfluid BCS-like state HI IH 113, the obtainment of the 
first samples of degenerate fermionic gases '28', '29^ , and 
some other recent developments i3Q^.^^,^,^. The- 
oretical studies have concentrated on the analysis of the 
BCS state [S^ [S^, and its excitation energies [39l |. 
as well as on comparisons between fermionic and bosonic 
properties j40ll4ll|. 

In the present work, we develop a dynamic theory of 
output coupling, for fermionic atoms initially confined 
in a magnetic trap. The outgoing atoms are considered 
as free particles. Our method can be easily generalized 
however to account for a gravitational field. We consider 
a one-dimensional model, with a spatially localized delta- 
type coupling between the atoms in the trap and the 
continuum of external modes. No external replenishment 
of the trap is considered, so that this model leads to a 
decay of the population in the trap and a non-stationary 
outgoing atomic beam. 

For an arbitrary confinement potential, we obtain gen- 
eral time-dependent expressions for the atomic operators 
corresponding to trapped and free atoms, the trapping- 
level populations, the spectral distribution and the first- 
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and second-order correlation functions of the atomic out- 
going beam. From these general expressions, we particu- 
larize the results for two important special cases: the infi- 
nite box and the harmonic oscillator. In order to identify 
more clearly the features associated to the fermionic na- 
ture of the trapped atoms, we compare our results with 
the corresponding properties for a bosonic beam. 

Of special interest is the infinite-time spectral distri- 
bution of the non-stationary outgoing beam. It differs 
markedly from the corresponding distribution previously 
calculated for single-level bosonic models, which leads to 
a single approximately Lorentzian peak Q^. For weak 
coupling, we find a distribution that reflects the structure 
of the discrete levels of the trap, exhibiting peaks that 
get higher and narrower as the energy increases. As the 
coupling gets stronger, the atom's spectral distribution 
is displaced towards a set of new energies that character- 
ize "dressed states" of the fermions in the trap. These 
dressed states result from the coupling between the un- 
perturbed fermionic states in the trap and the untrapped 
continuum. In the strong-coupling limit, there is a rein- 
forcement of the lowest-energy dressed mode, which con- 
tributes to the spectral distribution of the outgoing beam 
more strongly than the other modes. This effect is espe- 
cially pronounced for the infinite box, which indicates 
that the efficiency of the mode-reinforcement mechanism 
depends on the steepness of the confinement potential. In 
this case, a striking effect occurs: the fermionic beam be- 
comes quasi-monochromatic, in spite of the large number 
of energy levels populated in the initial trapped fermionic 
system. As expected, the fermionic second-order correla- 
tion functions exhibit the property of anti-correlation. 
Under these conditions, we predict therefore a quasi- 
monochromatic anti-bunched fermionic atomic beam. 

The paper is organized in the following way. In the 
next section we introduce the physical model and write 
down the basic Hamiltonian. The eigenvalue spectrum 
for this Hamiltonian is analyzed in Sec. IIIII where it is 
shown, for two forms of the trapping potential (infinite 
box and harmonic oscillator) that the one-dimensional 
model always exhibits a bound state, for any value of 
the coupling constant. In Sec. IIVI we analyze the non- 
Markovian behavior of the number of atoms in the trap. 
The spectral distribution of the outgoing atoms is ana- 
lyzed in Sec. General expressions for the field opera- 
tors and the correlation functions of the outgoing atoms 
are derived in Sec. I VII The corresponding numerical re- 
sults are presented in Sec. IVIII where comparisons are 
made between the bosonic and the fermionic cases. Some 
detailed calculations are referred to a set of two appen- 
dices. 



II. THE MODEL 

The physical model considers a single atomic species 
in a one-dimensional magnetic trap, with an external 
electromagnetic field inducing transitions between each 



trapped level and a continuum of non-trapped states. We 
ignore the effect of the magnetic field on the non-trapped 
state, assuming for instance that the trapped level corre- 
sponds to an electronic spin component -1-1, while the 
non-trapped state corresponds to the spin component 
zero. The total spin of the atom, nuclear plus electronic, 
is assumed to be a half-integer, so that the atom is a 
fermion and we neglect the small effects due to the nu- 
clear magnetic moment. We model this system by an 
effective Hamiltonian, with a bilinear coupling involving 
the field operators for the trapped and untrapped atoms. 
The one-particle eigenf unctions of the trapping potential 
are denoted by <^n(a;), the corresponding energy levels 
being given by ?iw„. The untrapped states are identi- 
fied by the center-of-mass wavefunction ip^^x), labeled 
by a continuous parameter ^, with energy Huj^. Thus, if 
we consider the untrapped atoms as free particles with 
mass M, we have hcu^ = (?i^)^/2M, and is the atomic 
momentum, while tp^i^) = exp(i^x)/'\/27r [normalized so 
that / dxip'^{x)ip^' (x) — — ^')]. If the atoms were un- 
der the action of a gravitational field, for instance, then 
^ would be an Airy-function index. The coupling is as- 
sumed to be spatially localized and is represented by a 
delta function. While this seems to be a most unphysical 
assumption, one can imagine a realization where a very 
tightly focused pair of Raman laser beams, in a weak 
trap, would induce the electronic transitions. As long as 
the beams' waists are much shorter than the de Broglie 
wavelength of the atoms, which could be the case for the 
lowest-energy states, the assumption of a delta-function 
coupling can be a good one. 

The effective Hamiltonian is written as 



H — Hx + Hp + He ■ 
In this expression, 

Ht = huJna'lan , 



Hp 



(1) 

(2) 
(3) 



where Ht corresponds to the trapped atoms, a„ is the an- 
nihilation operator for an atom in the trapping-potential 
eigenstate \ipn), Hp describes the untrapped atoms, and 

is the operator that annihilates an untrapped atom 
with wavefunction tp^^x). 

For fermionic atoms, the above operators obey the 
anti-commutation relations: 



{a„,SU = 0. 



(4) 
(5) 
(6) 



The coupling part of the Hamiltonian is given, for a 
general spatial-dependent coupling, by 



Hc^ih / dx\{x)i!\x)^{x) + H.c. 



(7) 
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where the field operators are given by 



(8) 



(9) 



In terms of the operators 6^ and a„, the interaction 
He may be written as 

He = ih^ j dign{^)h\an + YL.c, (10) 



where 



QniO = / dx \{x)lj)(^{x)(pn{x) . 



(11) 



In the special case of a delta-function coupling X{x) 
\6{x), Eq. (irU)) reduces to 



Hr 



ihX 



d^reiO)MO)bUn + il.c., (12) 



which is the interaction used throughout this paper. 
Without any loss of generality, the coupling constant A 
is taken to be real. 

This model may be considered as a multilevel extension 
of other systems considered before, which have a single 
level interacting with a continuum |42j |. 

For sufficiently strong coupling, bound states have 
been shown to appear in bilinear Hamiltonians involv- 
ing the interaction of a single mode of the electromag- 
netic field with a photon reservoir 141^ . The correspond- 
ing one-dimensional model, with free massive particles as 
the reservoir, was shown to exhibit a bound state for any 
value of the coupling constant ^2 ■ We may thus suspect 
that the above Hamiltonian also exhibits bound states. 
This is proven in the next Section, for two special cases 
of the trapping potential. 



III. DIAGONALIZATION OF THE 
HAMILTONIAN 

We generalize in this section the procedure adopted 
in ^3 for bosons in a single trapped level coupled to 
a reservoir of free massive particles. The existence of a 
multitude of bound levels in our case does not allow one 
to reach general conclusions concerning the existence of 
bound states for any trapping potential. We consider 
therefore two specific examples, the infinite-box and the 
harmonic oscillator, and show that the coupling given in 
Eq. H12|) leads to the existence of a single bound state, 
for any value of the coupling constant. 

We take the untrapped atoms as free massive particles 
(no gravitational field) , so that ■0fc (0) = 1/ V^n and the 



Hamiltonian of the system may be written as 

/ + 00 
dkhcokbih 
-oo 



ihX 



2tt 



+ 00 



^V3„(0)a„ / dkbl + U.' 



(13) 



where in this case 

As in |12| | , we introduce the even and odd operators ct 
and dk, given by 



1 



Cfc = — =[6fe + , 

dk = — =[-5fc + 6-fc] 
v2 



(14) 
(15) 



The operators dk are not coupled to the trap, so we 
can consider only the operators c^, and write 

H = y^hujndl^dn + / dkhukclck 
Jo 



^V<^„(0)a„ [ , 
Jo 



dk c\ + H.c. 



(16) 



In order to diagonalize this Hamiltonian, we apply 
Fano's procedure 0,^3 1 introducing the operators 



Ak = V'a„(A;)a„ -|- / ( 
„ Jo 



dk'-/{k,k')ck' , (17) 



so that 



iJ= / dkhn{k)AlAk, 



(18) 



where E^:^^ is the lower bound of H . 

Since the Hamiltonian given by Eqs. id, ©, ©, and 
()12|l is quadratic in the atomic operators, the same pro- 
cedure holds for fermionic and bosonic atoms. We start 
by calculating the commutator [Ak , H] , which yields two 
equivalent expressions, obtained by using either Eq. l(TH|) 
or Eq. d for H: 

[Ak,H]^hn{k)Ak 

r + co 

= hfl{k) Q!„(fc)a„ -|- / dk'j{k,k')ck> 
Jo 

= flUJnan{k)dn ~ ^ (y5* (0)a„(fc) 
^1 ^ n 

x: 

dk'hujk'^{k, k')ck' 

\ r+oo 

= V^„(0)a„/ dk'^ik,k'). (19) 
Jo 



X / dk' Ck 



ihX 
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From this equality we obtain the following equations: 
[n{k) - c^„]a„(fc) = -^M^) / dfc'7(fc, k') , (20) 



^^<,(0)a„-(fc). (21) 



We consider first the negative-energy solutions of these 
equations, which correspond to bound states. 



A. Bound states 

For bound states, we may set r2(fc) = — /i^, /i > 0, 
so that the bound-state energy is Eb = —^f^^, and let 
in this case a„(fc) ctfi.n, l{k,k') 7^(fc'), so that 
Eq. 117|l is replaced by 

POO 

= ^ a^,„a„ + / 7^(fc)cfc dk (22) 

and Eqs. ((201 and (|?T|l become: 

a^,„ = l^^dkj.ik), (23) 



iA 1 



^^:,(0)a^X- (24) 



n' 

Substituting Eq. ^ into Eq. (|^ . we obtain 
_ y^„(0) /-+°° dk 



(25) 



Multiplying the last equation by (fi'^iO) and summing 
over n, we have 



2 /--l-oo 



TT ^ fi'^ + LUn Jo fJ-'^ + t^fe 



dk 



(26) 



This equation immediately yields the eigenvalue equa- 
tion for /i: 



where 




27r 7o y + V 2% 

and we have used in Eq. that = hk'^/2M. 



(27) 

(28) 
(29) 



Replacing Eq. (^5)) into Eq. lf?7|l . we get the eigenvalue 
equation 



A^ 



2M 



(30) 



We can see that this equation has one and only one 
solution if F{y) is finite when y — 0, and F{y) — > 
when y — !■ oo. This will be shown to be the case for the 
two special cases considered in this paper. One should 
note however that the form of this eigenvalue equation is 
highly dependent on the dimensionality of the problem. 
This dependence is quite apparent in the expression for 
the function /(y), where the divergence for y = disap- 
pears if one replaces in Eq. H29|l dk by d^fc (adding up 
a cutoff to the upper integration limit, so that the inte- 
gral remains finite). This would imply the replacement 
of fi on the right-hand side of Eq. H30|) by a function of 
fi that would not go to zero when fj, —f 0, and there- 
fore the bound state would appear only for a sufficiently 
strong coupling. For bosons at zero temperature, the de- 
pendence on the dimensionality of the bound state of the 
corresponding Hamiltonian (with just one bound-level) 
was explicitly demonstrated in Ref. [T2| . 

The functions a^_„ and 7^ (A;) may be obtained in the 
following way. We impose the condition {A^, Aj^^} = 1 
(for bosons we would replace the anti-commutator by a 
commutator, with the same results at the end), obtaining 



+00 



dk\-f,{k)\' = 1 . 



(31) 



Replacing Eqs. iPUl and (^51) into this equation, we 
obtain, except for an irrelevant overall phase factor that 
can be absorbed into the definition of the states (pn{0)' 

^;(0)a,,„ = o^yL,. , (32) 



where F'(/i^) is the derivative of F{y), defined by 
Eq. |(2H|), evaluated at y = ^^: 



(33) 



Taking Eq. ^ into Eqs. ^ and (gSJ), we get finally 

_ lpn{0)^i/i^i'^ +UJn) 



(34) 



(35) 



We discuss now the solutions of Eq. (|30|l for two im- 
portant special cases of trapping potential: the infinite 
box and the harmonic oscillator. We show that in both 
cases there is one and only one bound state, for any non- 
vanishing value of the coupling constant. This implies 
that Eq. I|18() becomes 



H= dkhujkA^k)A{k) -hfi^AlA^, 



(36) 
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where cok = hk'^/2M, —hfj,"^ is the energy of the bound 
state, and is the corresponding annihilation operator, 
given by Eq. . 



1. Infinite box 
In this case, we have for the trapped particles. 



2 2 
-n — LUifl , 



(37) 



2ML2 

the corresponding eigenstates being given, for odd n, by: 



— cos (fc„r) , 



(38) 



where L is the length of the box. 

For even n, the cosine function is replaced by the sine 
function, which vanishes for x — 0. The interaction in 
Eq. l(T^ does not couple these states to the outgoing 
beam, and they do not contribute to the sum defining 
F{^'^) (this is a consequence of the localized nature of 
the symmetric coupling). Thus, only odd n's (even wave- 
functions) contribute to F{^'^) , which may be written as: 



F 



(box) 



2 

uJiL 



E 

n (odd) 



1 



(39) 



From Eq. HA9|) of Appendix fXl we have 



F 



(box) 



tanh {ji[il1yfuj\ 



(40) 



Therefore, in this case F{iJ,'^) goes to a finite value when 
/i — > 0, and vanishes when fj, oo. It is clear then 
that Eq. H30|l has a unique solution. In terms of the 
adimensional coupling constant S defined by 



S — Xtt'^ /Lull , 



(41) 



the weak-coupling limit corresponds to (5 ^ 1, so that 
<C oJi, and the hyperbolic tangent may be approxi- 
mated by its value close to the origin, thus yielding 



Ef 



(42) 



In the strong-coupling limit (5 3> 1, we get /i^ 3> wi, 
so that the hyperbolic tangent may be approximated by 
one, and 



Ef 



2tt' 



:huji 



(43) 



2. Harmonic oscillator 



For the harmonic trap, we have 

flLOn = hLUo{'n -1-1/2) 



(44) 



and 

where Hn{x) is the Hermite polynomial of order n and 
d — \/h/ mujQ is the width of the ground state. 

Wave functions corresponding to odd values of n do 
not contribute to the sum defining F{p'^), which becomes 
now: 



|^2m(0)P 



2^0 TO + 1/4 + ifi^/2uJo) 

= ^ , r(i/4 + A^V2^o) 

2uJod T{3/4 + fiy2LUo) ' ^ ' 

where r(a;) is the Gamma function This result is 

proven in Appendix IbI 

Replacing Eq. H4t)|) into Eq. (|30|l . we obtain the final 
expression for the eigenvalue equation 



2y/2uJod , TTfT. — 7 = M ■ 



r(3/4 + /x2/2c^o) 



where 5' is defined by 



5' = \/2ujQd. 



(47) 



(48) 



It is easy to verify that Eq. H47|l has one and only one 
solution /i > for any 5' . 

In the weak-coupling limit 5' <C 1, we may neglect 
the contribution of /i in the argument of the Gamma 
functions, thus getting 



Eb = -h^? = 



r(i/4) 



n 2 



_r(3/4) 

For strong coupling, we use the following identity ^ 



(49) 



r(6 + z) 



l + 0(z-i). 



(50) 



Identifying z /2ujq, a — > 1/4 and b 3/4, we 

obtain for 6' ^ 1: 



Eb = -hfi^ = . 



(51) 



B. Positive-Energy Solutions 

From Eqs. and we can write for Vl{k) = ujk> 



/TT UJk — UJn Jo 

iX 



dk'-f{k,k'), (52) 



UJk — UJk' 



+ Z{k)S{uJk - UJk') 



(53) 
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where P stands for the principal part, and we have as- 
sumed for the moment, in getting Eq. (|52|l . that a;^. ^ a;„, 
for any n. In these equations Z{k) is a function to be de- 
termined. Inserting Eq. 152|) into Eq. (|53|) . we obtain the 
expression for Z{k): 



Z{k) = 



fik 



(54) 



where F{y) was defined in Eq. 1281) . and we have used 
that S{u>k -ujk') = {M/h\k\)6{k - k'). 

Using Eqs. ((T7|l . (|^ . (|^ . and lf5j|) . and imposing the 
condition 



{i(fc),it(fc')} = <5(fc-fc'), 



we obtain 
Therefore, 

and 



^<(0)a„(fc) = 



-iihk/M) 



nk 



M Xy^TT^ + Z^{k) 



(55) 



(56) 



Wfc - UJk' 



Z{k)5{ujk - UJk') 
(57) 



a„(fc) 



A (p„(0) Z(fc) 



/TT Wfe - a;„ ^7r2 -)- Z2(fc) 



(58) 



From Eqs. H28|l and (|54|l . it is easy to check that a„(fc), 
given by Eq. H58|l. remains finite when Wfe ^ u;„. This 
allows one to remove the restriction Wk ujn, used to get 
Eq. (|52|l . and adopt Eq. (|58|l as the expression for Q:„(fc) 
for all values of k. 



IV. POPULATION IN THE TRAP: 
NON-MARKOVIAN BEHAVIOR 

One of the consequences of the existence of the bound 
state is the failure of the Born-Markov approximation for 
this problem. A related consequence is that a fraction 
of the atoms remains in the cavity, even in the infinite- 
time limit t ^ oo. This can be seen by writing down 
the decomposition of each cavity mode in terms of the 
eigenmodes of the Hamiltonian: 

/"OO 

a„(t)=/ <(fc)ifce-*'"''='*/2*^dfc + a* e'^'*4, (59) 



and replacing the operators Ak and by their expres- 
sions in terms of the operators a„ (0) and Ck (0) . One gets 
then: 

a„(i) = V[ /"<(fc)a„,(fc)e-^'"='*/2M^^ 

n' L^O 



dk' 



rffc<(fc)7(fc,fc')e-^'"'"''/'*^ 



/\ ifj, t 



Cfc'(O) 



(60) 



This expression exhibits explicitly the coupling be- 
tween the trap modes, which is induced by the coupling 
with the external modes. If initially only the trap modes 
are populated, and if one is interested only in normal- 
ordered correlation functions, the contribution of the op- 
erators Cfe(O) may be ignored. This will be always the 
case in the present paper. 

If at time t = only the cavity mode n is populated, 
the fraction of atoms left in the same mode at a later 
time t is given by 

7t 



{aiit)anit)) 
(ai(0)a„(0)) 



|a„(fc)pe-''"='*/2M^fc-K|a, 



.(61) 



The integral vanishes in the infinite-time limit, since 
Q!„(fc) remains finite for all values of k, and therefore 
in this limit the fraction of atoms left in the cavity is 
This result is easy to understand: in order to get 
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the residual population, one must multiply the fraction 
of atoms in the initial mode that are in the bound mode 
A^, given by ja^i.nP, by the fraction of the cavity mode 
fln present in A^, which is also given by |q;^,„P. Fur- 
thermore, the time-dependent population exhibits oscil- 
lations, resulting from the beating between the integral 
and the discrete contribution in Eq. (|61|l . 

A similar behavior holds if initially more than one 
bound mode is populated, as it is the case for trapped 
fermions at zero temperature. The residual population 
of level n is then given by 

{alan)i^) - |a^,„P^ (al,(0)a„.(0)) , (62) 



where now J2n' l'^M,n'P(^I'(0)'Jri'(0)) is the fraction of 
the initial population that is in the bound mode , and 
|q^m>"P is fraction of the cavity mode a„ present in 
the bound mode. 

An expression for the total residual population N{oo) 
inside the trap may be obtained from Eqs. (|62() and 134|l : 

N{^) = ^(ata„)((^)=5]|«^^„f (4,(0)a„,(0)) 

n n' 



(63) 



where F'{iJ?) is, as before, the derivative of F(y), given 
by Eq. evaluated a,t y = fi"^. 

For TV bosons at zero temperature, only the term 
with n' — contributes to the above sum, and 
(a]^, (0)a„'(0)) — N. On the other hand, for fermions at 
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zero temperature, the initial population is the same for 
all levels (one atom for each level, since all the trapped 
atoms have the same spin), up to the last occupied one 
(Fermi surface). If the number of atoms is much larger 
than one, then one may approximate the sum in the 
above expression by one with an infinite number of terms. 
The resulting number is an upper bound for the residual 
population inside the trap, which is actually achieved 
when N oo: 



:(oo) 



(64) 



In the weak-coupling limit (S <C 1), an approximate 
expression for H63() may be obtained, for the 1-D box and 
the 1-D harmonic oscillator, by using the results obtained 
before for the function -F(/i^) and for the bound-state 
energy Eb ~ —h^P. We get thus, for the 1-D box. 



N{oo) 



96 ^ 

i(odd) 



'at(0)a„(0)) 



(65) 



and for the harmonic oscillator 



iV(oo) = 64:V^S'^ 

oo 

xE 



r(i/4) 



r(3/4) 

(2n)! 



n=0 



22"(n!)2(n+ 1/4)^ 



(4„(0)a2„(0)). (66) 



From these results, we can see that the residual popu- 
lation inside the trap is very small in the weak-coupling 
limit, being proportional, both for the box and the har- 
monic oscillator, to the eighth power of the corresponding 
dimensionless coupling constant. 

In the strong-coupling limit {6 3> 1), we get both for 
the box and the harmonic oscillator that the upper bound 
for the population inside the trap is -/Vniax(oo) — 1/4. 
This result, which is actually achieved when the number 
of atoms is much larger than one, shows that a substantial 
fraction of the atoms remains in the trap in the infinite- 
time limit. This is a direct consequence of the existence 
of a bound state of the total Hamiltonian. 

We proceed now to the calculation of the spectral dis- 
tribution of the outgoing atomic beam. 



V. SPECTRAL DISTRIBUTION OF THE 
OUTGOING BEAM 

A time-dependent spectral distribution for the outgo- 
ing fcrmionic beam can be obtained from the expression 
of the free-atom operators in terms of the operators that 
diagonalize the Hamiltonian: 

Ck{t)^ / dfc'7*(fc',A:)e-'-^-'*ife,(0)-f7;(^)e^^'*4(0)- 
Jo 

(67) 



From Eqs. |(T7|), IHl), and (gZI), we get, ignoring the 
contribution of the operators Cfc(O) (since initially the 
outside modes are empty, and only normal-ordered cor- 
relation hmctions are considered): 



Cfe(i) 



E 



dk'j*{k',k)an{k')e 



— iuj^i t 



7u(fc)ap 



a„(0). 



(68) 



The time-dependent spectral distribution is given by 
(blit)bk{t)) . This quantity can be expressed in terms of 
Ck{t), using that from Eqs. H14|) and (|15|l . 

1 

1 



b. 



--{ck - dk) , 
= (cfe + dk) , 



(69) 



and that dk does not couple with the trapped-atoms oper- 
ators, so that it can be ignored when calculating normal- 
ordered correlation functions [of course, its presence in 
Eq. H69|l is important to get the correct commutation 
relations for the operators bk and Ck]- One gets then 
that the time-dependent spectral distribution is given by 
{ciit)ckit))/2. 

The integral in Eq. H68|) can be calculated by using 
Eqs. (|57|l and (|58|l . For finite times, one has to con- 
sider the contributions from the complex poles of the in- 
tegrand, which give rise to exponentially decaying terms. 
These contributions can be handled numerically. An ex- 
ample will be given in Section IVIII 

An analytic expression can be obtained in the infinite- 
time limit. Since the contributions from the complex 
poles of the integrand in Eq. I|68|l give rise to exponen- 
tially decaying terms, they will be negligible in this limit, 
so the relevant contributions come from the principal part 
and the delta function in Eq. lf57|) . One gets then: 



Cfe(i) 



E 



-A(^„(0)/V¥ 



LUk-LUn X^{M/hk)F{-UJk)+i 



an(0) 



(70) 



We will show in the following that the first term on 
the right-hand side of this equation leads, in the weak- 
and strong- coupling limits, to narrow peaks, centered 
around the unperturbed energies of the trapped atoms in 
the weak-coupling case, and around dressed energies of 
the coupled atoms in the strong-coupling limit. The last 
term on the right-hand side of this equation is the bound- 
state contribution. It can be shown to be much smaller 
than the remaining terms in Eq. (|70|l in the regions of 
the spectral distribution close to the peaks, so it will be 
neglected from now on. One gets then, using Eq. 169() 
and neglecting the contribution from dk'- 



E 



-A^„(0)/V2^ 



UJk-UJn X'^iM/h\k\)Fi-LOk)+i 
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X a„(0). 



(71) 



In view of Eq. (|69|l . this equation yields bk for both signs 
of fc. 

From this expression, and assuming that the initial 
state is diagonal in the number representation, one gets 
the outgoing beam spectral distribution in the infinite- 
time limit: 

X (at (0)a„(0)) . (72) 
From Eqs. H54() and (|58ll . it is easy to see that 

(&l6fe)(oo) - i ^ |a„(fc)p(at (0)a„(0)) , (73) 

n 

which shows that the contribution to the spectral distri- 
bution from trap level n is proportional to the probability 
|Q;n(fc)P that the atom in this level is in the eigenmode 
of the total Hamiltonian with energy cuk . 

We consider now the specialization of Eq. 1)72(1 to the 
infinite-box and harmonic oscillator potentials. 



Infinite Box 



From Eq. we get: 



j.(box)( ^ tan(fcL/2) 



(74) 



Taking this result, plus Eqs. 1(23, (EHIl, and into 
Eq. H71|) . we get, for large times: 



^ (fcL)2 

=i(odd) ^ ' 



5./L/^{kLf 



.(0), (75) 



where 



G(fc,i) 



-cos(|fcL/2|)e-^'*'''*/2M 



i{kLY cos (|fcL/2|) (52/4) sin {\kL/2\) ' 

(76) 

One should note that the singularities in the sum are 
canceled out by the numerator of G{k, t). 

We discuss now the behavior of these expressions in 
two limiting cases, corresponding to weak [5 <^ 1) and 
strong {S ^ 1) coupling. For weak coupling (5 <C 1), 
the term proportional to cos(|fcL/2|) dominates in the 
denominator of G{k,t), which exhibits sharp resonances 
close to values of k that correspond to the bound states 
of the infinite-box potential: kL = nir, n odd - these 
are the zeros of cos(|fcL/2|). We may thus approximate 
the expression in Eq. (|75l) by setting, around each peak, 
kL = nn + Pn, \0n\ ^ 1, keeping only the lowest-order 
terms in the expansions of the trigonometric functions 
in Eq. I|7t)|) . and neglecting small corrections in /3„ for 



the other fc-dependent contributions. Neglecting these 
corrections means that deviations from the Lorentzian 
shape will be ignored here. 

One gets then, from Eqs. H75|) and (|76|) . that, asymp- 
totically in time. 



where 



^(„) _ i(5/L7^e-*'"''='*/2^V2™ 



[kL - riTr) + 2i{S/2n7ry 



a„(0). 



(77) 



(78) 



The infinite-time spectral distribution is then given by 
the sum of the contributions from all peaks: 

mioo)/L= E W ^kL^nn/+i5/V2nn)^ 
n(od.d) 



X (at(0)a„(0)) 



(79) 



Under the above approximations, each term in the sum 
is a Lorentzian with width (AfcL)p-y^jj]y[ — 2{S/mr)'^, 
and height increasing with the square of n. Therefore, 
the peaks in the spectral distribution become higher and 
narrower as n increases, in such a way that the area 
under each peak is proportional to the state population 
(fit (0)a„(0)) corresponding to the same value of n. One 
should note that the contribution for each peak stems, in 
this case, from a single bound state. The total number of 
peaks is equal to the number of initially populated states. 

For strong coupling, 5 3> 1, the term proportional to 
sin(|fc^/2|) dominates in the denominator of Eq. (|76|l . im- 
plying that the peaks are shifted towards the values of 
k corresponding to the zeroes of sin(|fcL/2|). These val- 
ues, given by kL = 2m7r, m integer, are precisely mid- 
way between the weak-coupling peaks, and are associated 
to dressed energies of the system, which originate from 
the interaction between the discrete states through the 
continuum, as shown in Eq. ((6()|l . Indeed, it is easy to 
see that, in this case, the sum over n must be kept in 
the approximation analog to Eq. (|78|l of the expression 
in Eq. (|75(l around kL = 2m7r, m integer. Therefore, 
several trap states contribute now to each peak in the 
spectral distribution. 

The complete infinite-time spectral distribution is 
given, as before, by the sum of the contributions for all 
peaks: 

^ (4V2m7r/5)4 + {kL - 2m7r)2 



E 

n=i(odd) 



(4to2 — n2^2 



'4(0)a„(0)).(80) 



This expression differs remarkably from the one in 
Eq. (|79|l . There is now an infinite number of peaks. 
Each peak, with width (AfcL)p^jj]yj = 32(m7r/(5)2, is 
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now fed by all the populations in the trap. For the peak 
at kL — 2m7r, the strongest contributions come from 
the two populations with n = 2m ± 1 (first neighbors of 
kL = 2rmr), with heights proportional to 1/(4to + 1)^. 
The contributions from the second neighbors is propor- 
tional to l/9(4r7i + 3)^, which is less than nine times 
smaller. The contributions from states that are farther 
away decrease as l/n"*. This implies that the height of 
the peaks decreases as 1/m^, while its width increases 
as m^, and that the contribution for each peak comes 
mainly from the populations of the first neighboring trap 
states. These two states are the main components of the 
dressed state that contributes to this peak. 



B. Harmonic Oscillator 

Setting fM^ = —LOk in Eq. (|^ . we get: 

^(ho)(_^^^)^^£lVi^(MM 



2ujod r[3/4- {kd/2y 



(81) 



Inserting this result, plus Eq. 144|) into Eq. H71|) . we 



get: 



where 



Jiv.t) 



S'd ^ ip2rr,{0)\kd/2\J{kd,t) 

(to + 1/4) 



^ {kd/2y 

m— 



a2m(0), (82) 



r"i(l/4-2;V4)2)e-*^'='*/2A'^ 



(5'2r-i(3/4 - yV4) + i|?//2|r-i(l/4 - ' 

(83) 

The spectral distribution of the outgoing beam is given 
by {blbk){oo). From Eq. JHU, we can see that, for 6' < 1, 
the peaks of the spectral distribution should be close to 
the zeroes of r"i[l/4-(fcd/2)2], while for ^' > 1 they ap- 
proach the zeroes of r~^[3/4 — {kd/2)'^], which are, as in 
the infinite-box case, midway between the weak-coupling 
peaks. These zeroes correspond to the dressed energies 
of the harmonic trap, under strong-coupling conditions. 

We obtain now approximate expressions for the spec- 
tral distribution in the weak- and strong-coupling re- 
gions, by expanding the Gamma functions about the val- 
ues of k that correspond to the spectral peaks. A useful 
equality for this purpose is 



T-\x) = r(i - x) 



sin (ttx) 



(84) 



For weak coupling, the peaks are around {kd/2)'^ — 
m+ 1/4, TO integer. We get then, by expanding the func- 
tions of k in Eq. (|82|l around these values, using Eq. (|45|l . 
and approximating the spectral distribution by the sum 
of the contributions from all the peaks: 



{blbk){^)/d 



(2m)!/22"(TO!)2 



X (aL(0)a2™(0)) , 



- 1/4] 



(85) 



where the linewidth Tm is given by 
S'^ r(TO-f i) 



r — 



(86) 



For TO ^ 1, we find, using Stirling's approximation, 
that the linewidth of the peak of order to is given by 



r = 

m. 



(87) 



while the corresponding height is ttv^^^ /2n^S'^ . 

In the strong-coupling regime the peaks are around the 
values (fcd/2)2 = to -I- 3/4, so that one gets, approximat- 
ing the contribution around each peak: 



{bibk){^)/d = 



(to + |)r2(TO-K 

2J'27r7/2 (TO-n + i)2 

m— n— ^ ^ ^ 

(2n)!/22"(n!TO!)2 



[(M/2)2-TO-|]2-f f', 

x(aLa2n(0)) , 
where the linewidth is now given by 



r' — 



(89) 



One should note that, in the strong-coupling limit, all 
trapped-level populations contribute to each resonance, 
as opposed to the weak-coupling limit, when each res- 
onance is associated with a single trapping level. The 
same phenomenon occurred in the infinite-box potential. 

In this case, for to 3> 1, we have for the linewidth of 
the peak of order to. 



TO 

6^ 



(90) 



while the corresponding height is proportional to 



VI. FIELD OPERATORS AND CORRELATION 
FUNCTIONS 

A. Field operators 

The time-dependent field operators for the outgoing 
atoms are given by: 



^ix,t) = I dk^bkit). 
1 V 27r 



(91) 



From Eqs. H68() and (|69|l . this can be written in the fol- 
lowing form, if we ignore the vacuum terms proportional 
to Cfc(O) and dfc(O): 



*(x, t) = ^ N{n, X, t)a„(0) , 



(92) 



10 



where 

N{n,x,t) 



X e 



+°° dk_ 
\^ 



dk'j*{k\k)an{k') 



Akx 



(93) 



We may call this the source contribution to the field 
operators. One should note that iV(n, x, 0) = 0, as ex- 
pected (no contribution, at i = 0, of the trap modes to 
the field operators corresponding to the outgoing atoms). 
Indeed, from Eqs. H59|) and (|67l) . {cfe, } = implies that 



dfc'7*(fc', k)an{k') + 7:(fc)aM,n = . 



(94) 



From Eqs. 1)92(1 and (|93|l . we can see that the field oper- 
ator corresponding to the outgoing atoms is given by the 
sum of two contributions, besides the terms proportional 
to Cfc(O) and dfe(O), 

i>{x, t) = ^(bound) ^) ^ ^(run) ^) ^ ^g^^ 
where the bound-state contribution is given by 

^^;[k)a,,^e'-e^^'' 

(96) 



^(bound)(^^^^ ^Y,l 

X a„(0), 
and the running-wave part is 

'■+°° dk 



(run) 



dk'j*{k',k)an{k') 



X e~'"'='*cosfca:a„(0) . 



(97) 



The bound-state contribution is readily calculated, by 
using Eqs. ^ and 



^(bound)^^^^) 



E 



F(^2) _ ^2^/(^2) 
</3„(0) 



-an(0) , 



(98) 



which exhibits a spatial dependence that decays expo- 
nentially, with a decay constant given by y/2M\EB\/Ti, 
where Eb — —hji^ is the bound-state energy. 

For the running- wave part, simple results can be ob- 
tained by replacing directly into Eq. (|91|l the asymptotic 
results obtained for the operators bk in the weak- and 
strong-coupling limit, for the special cases of the infinite 
box and the harmonic oscillator. We restrict ourselves 
here to the infinite-box case, in the weak-coupling limit, 
since the results for the strong-coupling limit and the 
harmonic oscillator are quite similar. 

Taking Eqs. G?)) and lO into Eq. (jHIl), one gets, in 
the weak-coupling limit: 



#(run)(^^,)^^ ^ iv('-H(n,a;,i)a„(0) 



.(odd) 



(99) 




FIG. 1: Magnitude iVi(r) = \N^^^^\i,x,t)\ of the outgoing 
wave packet corresponding to the infinite-box level n = 1, at 
the position x — 20L, where L is the box width, as a function 
of the renormalized time r — Tit/2ML^ . 



where 



iv(™'i)(n,x,r) 



-t-oo 



dy e cos{yx/L) 



n [{y -rnr) + 2i{S/2mr)^] ' 

(100) 

and T is the renormalized time r — M/2ML'^. 

Equations and (|100|l show that the outgoing field 
corresponds to a train of wave packets, centered in mo- 
mentum space around the momenta corresponding to 
the trap eigenenergies. The expression in E g. (I100|l is 
closely related to the Moshinsky function 03 > which 
yields the time-dependent behavior of a wave packet ini- 
tially confined in a half-space: 'ip{x, 0) — d{—x) exp{ikx), 
3m fc < 0, where 6{x) is the Heaviside function [9{x) = 
for a; < 0, 6{x) = 1 for x > 0]: 



Mix,k,t) = — 



+ 00 



dhi- 



piKX — hK, t/2ri 



(101) 



This function exhibits an oscillatory behavior as a func- 
tion of time, the so-called "diffraction in time" ef- 
fect IItIEL This effect, which has been experimentally 
observed [43|, appears when a shutter placed at a; = is 
opened, letting the initial wave packet, confined to the 
X < region, evolve. Here however the integration is 
from to infinity, the difference stemming from the fact 
that in our case the outgoing field emerges from a; = 
and propagates in both directions. In spite of this differ- 
ence, we also get here transient effects that can be de- 
scribed in terms of a diffraction in time. This is clearly 
shown in Fig.^ which displays the plot, as a function of 
the renormalized time r, of the magnitude of the expres- 
sion in Eq. ((TUni) for n = 1, a; = 20L, and S^/2tt^ = 0.2. 
In this figure, the small-amplitude fringes close to the ori- 
gin are due to the interference between the wave packet 
propagating towards the positive direction with the tail 
of the packet that propagates in the negative direction. 
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This tail is present even at t = 0, contrary to what hap- 
pens with the Moshinsky function, which vanishes ex- 
actly for a; > at the initial time. Here, when t — 0, the 
sum of the tails of all the wave packets exactly cancels out 
the bound-state contribution, leading to the vanishing of 
the source contribution at the initial time. 



B. Correlation functions 

From Eqs. (|92|l and H93|l . we can calculate the normal- 
ordered correlation functions of any order for the system 
in question. Other orderings could also be considered, 
by keeping the terms dependent on 6^(0) in the above 
expressions. 

Thus, for an initial state diagonal in the number rep- 
resentation, the first-order coherence is given by 

Gi{x,x',t) = {^''{x,t)^ix',t)) 

= ^ iV* (n, X, t)N{n, x' , t)(at a„) (0) . (102) 

n 

For X = x' , this becomes the beam density I{x, t). The 
normalized first order correlations functions is defined as 



gi{x,x\t) 



Gi{x, x' , t) 



(103) 



^I{x,t)I{x',t) ■ 

The normal-ordered second-order correlation function, 

G2{x, x, t) = {¥{x',t)¥{x, i)*(x, t)*(x', t)) , (104) 

may be written, for fcrmionic atoms, as '5C|: 

G2{x,x',t) = I{x,t)I{x\t) - |Gi(x,x',t)p. (105) 

The normalized second-order correlation function is 
defined as: 

(106) 

The minus sign in the above expression accounts for the 
anti-bunching property of fcrmionic beams. This second- 
order correlation function is analyzed numerically, for the 
special cases considered in this paper, in Sect.ltnn 

For bosonic atoms at zero temperature, the second- 
order correlation does not depend on the position, while 
for a thermal distribution it exhibits the bunching effect 

M. 



VII. NUMERICAL RESULTS 

We compute now the behavior of the output atomic 
beam, for zero temperature. We study the spectral dis- 
tribution of both a fcrmionic and a bosonic output beam, 
and also the second-order correlation function of the out- 
going fcrmionic beam. 

Our main numerical results are shown in Figs.|21to0 



o 

X 



A 



V 



o 

X 



A, 



V 



A 



3 

2 

1 • 

■ 
10 
8 

6 - 
4 
2 

5 
4 
3 



2 

V 1. 



(a) 





(b) 




(c) 



10 
kLV2 



15 



— I 
20 



FIG. 2: Normalized spectral distribution {h\hh)/L of the out- 
going fermionic beam, for an infinite-box potential of length 
L, for different coupling strengths: (a) 5 = 0.1; (b) 5 — 10; 
(c) 5 = 100. Initially, there are 21 atoms in the trap, corre- 
sponding to 11 coupled levels. The peaks in (a) are centered 
around the wave numbers corresponding to the energy levels 
inside the trap. As the coupling strength increases [Fig. (b)], 
the peaks are displaced towards the situation displayed in 
Fig. (c), where the resonances are related to "dressed states" 
of the system trap plus environment. While for weak coupling 
the height of the peaks increases with the energy, the opposite 
happens in the strong-coupling limit, when the lowest-energy 
peak is higher than the others. In these pictures only positive 
wave numbers are shown, since the complete graphic is sym- 
metric with respect to kL/2 = 0. Also, the totality of eleven 
peaks in the weak-coupling case is not displayed in Fig. (a). 



A. Spectral distribution of the output beam 



We consider first the spectral distribution 
(6j.(oo)6/c(oo)) of the output beam in the long-time 
limit. In the following, we take the initial number of par- 
ticles inside the trap to be iV = 21, so that the sums over 
the contributions of the trap levels will involve eleven 
terms (since only even states couple to the continuum). 
The actual shape of the spectral distribution depends 
strongly on the value of the adimensional parameters 
S and S', defined respectively by Eqs. and 
The spectral distribution for fermions trapped by an 
infinite box is displayed in Fig. [3 while the distribution 
corresponding to a harmonic potential is displayed in 
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FIG. 3: Normalized spectral distribution {b\bk)/d of the out- 
going fermionic beam, for a harmonic potential, where d is 
the ground-state width, for different coupling strengths: (a) 
S' = 0.1; (b) S' = 1; (c) 5' = 0.1. The behavior is similar to 
the one in the previous figure. Fig. (a) displays eleven peaks, 
the first one being barely visible and the last one correspond- 
ing to the Fermi surface for the = 21 trapped atoms. In the 
strong-coupling regime, one gets instead an infinite number of 
peaks with decreasing amplitude. The relative importance of 
the first peak is less pronounced here than in the infinite-box 
case. 
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FIG. 4: Normalized spectral distribution {h\.hh)/L of the out- 
going bosonic beam, for an infinite-box potential of length L, 
for different coupling strengths: (a) 5 — 0.1; (b) S — 10; (c) 
S — 100. For weak coupling - Fig. (a) - there is only one peak, 
which corresponds to the initially occupied trap level. As the 
coupling strength increases, new peaks appear, although only 
the first one remains relatively important. As before, this is 
an effect of the coupling of the trap levels through the con- 
tinuum. In this figure only the positive wave numbers are 
shown, as in Fig. 2. 



Fig.El 

We can see that for small values of 5 and 6', the peaks 
in the spectral distribution can be simply interpreted as 
resonances associated with the unperturbed trap levels 
corresponding to even eigenstates, consistently with the 
previous analysis. One should notice that the height of 
the peaks grows with k (the widths of the peaks cannot 
be seen, within the scale of the figures). In fact, in the 
weak-coupling limit the approximate expressions given 
by Eqs. H79|l and (|85l) lead to excellent fits to the spectra 
exhibited in Figs. [3 and O 

For larger values of 6 and S' , the peaks are shifted to- 
wards the strong-coupling dressed energies, as discussed 
in Sections lV Al and IV Bl Another important feature, also 
discussed in those Sections, is the relative amplitude of 
the peaks. We see that, in the strong-coupling limit, the 
first peak is the highest one, the amplitude of the succes- 
sive peaks decreasing now with the energy. This feature 
is especially pronounced for the infinite box. In this case. 



we have the remarkable feature that the outgoing beam 
has an accentuated monochromatic character, in spite of 
the multitude of bound levels of the infinite box. Again, 
there is excellent agreement of the strong-coupling plots 
with the approximate expressions given by Eqs. H8U|) and 
(Ell- 
in Figs. 21 and [S] we study the behavior of the corre- 
sponding bosonic systems, initially in the ground state. 
For small values of 6 and 6', the peak is related to the 
only occupied level of the system. For stronger couplings, 
new peaks appear, around energy levels corresponding 
to the excited trapped states, which get populated as a 
consequence of the interaction between the ground state 
and the continuum. We note that this effect is absent in 
models that treat the bosonic system as a single trapped 
level. 

The results shown so far correspond to the spectral 
distribution of the outgoing atoms in the infinite-time 
limit. The behavior of the time-dependent spectral dis- 
tribution for finite times is very simple in the case of 
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FIG. 5: Normalized spectral distribution {b\bk)/d of the out- 
going bosonic beam, for a harmonic potential, where d is 
the ground-state width, for different coupling strengths: (a) 
5' = 0.1; (b) 6' = 10; (c) 5' = 100. Fig. (c) clearly dis- 
plays both the energy displacement and the emergence of new 
peaks, in the strong-coupling case. 
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FIG. 6: Normalized spectral distribution {h\hh)/L of the out- 
going fermionic beam, for an infinite-box potential of length 
L, in the strong-coupling regime [5 = 100), for the renormal- 
ized time r = nt/2ML^ equal to: (a) r = 0.5; (b) t = 2; 
(c) r = 5. For r = 10 one recovers the infinite-time spectral 
distribution displayed in Fig. 2(c). 



weak coupling: then, each level behaves independently 
of the others, so that as the population of each trap level 
decays, as discussed in Section HVI the population of the 
corresponding free-space mode increases, finally getting 
to the distributions displayed in Figs. Inland 13 

For strong coupling, the dynamics of the spectral dis- 
tribution is more involved, since now the different trap 
levels interact with each other through the continuum, 
and there is a strong probability, for finite times, that 
the atoms are reabsorbed into the trap (this phenomenon 
shows up in the oscillatory behavior of the atomic popu- 
lation in the trap) . Figure 6 exhibits the approach to the 
infinite-time limit of the time-dependent spectral distri- 
bution, for the special case of strong coupling, and for 
an infinite-box potential. One should note that the rel- 
ative heights of the spectral peaks are time-dependent, 
their precise relationship being related to the intricate 
transient behavior of the system. 



B. Second-order correlation for fermions 

Figure [7| displays the behavior of the second-order cor- 
relation function as a function of the distance, for the 



fermionic case . The anti-correlation of fermions is clearly 
exhibited in all cases. For strong coupling, this correla- 
tion function goes very fast to one. 



VIII. CONCLUSIONS 

We have shown that the dynamics of trapped fermions 
with an output coupling may exhibit very interesting fea- 
tures. The main source of these features, which is also 
the most challenging aspect of this problem, as compared 
to the corresponding situation for bosons at zero temper- 
ature, is the multitude of energy levels of the trap that 
are necessarily populated, due to the Pauli principle. 

In view of the complexity of the problem, our strategy 
in this paper was to deal with a model simple enough 
so that it was possible to obtain some analytic handle 
on it, yet sufficiently rich to demonstrate interesting fea- 
tures of this system. The assumption of a delta-type 
coupling indeed greatly simplified the solution of the 
problem, while still keeping the main feature of leading 
to a coupling between the trap eigenstates mediated by 
the continuum. This coupling has remarkable effects in 
the strong-coupling limit, leading to clear signatures of 



14 



of 



1.0- 

0.8- 

0.6- 

. 0.4- 

0.2- 

0.0- 

1.0- 

0.8- 

0.6- 

0.4- 
i 

0.2- 
0.0- 



0.0 



R 



0.2 



of 



(b) !3> 







0.0 



0,4 



2 



4 



6 



8 



10 



R 



FIG. 7: Normalized second-order correlation function of the 
output fermionic beam as a function of the dimensionless 
position R: (a) Harmonic oscillator (7? = x/d, d being the 
ground-state width); (b) Infinite-box potential [R — x/L, L 
being the box length). The dimensionless time r, defined 
as r = ht/2ML^ for the box and as r = 2ht/Md'^ for the 
harmonic oscillator, is taken equal to 10. The full-line curves 
(shown in detail in the insets) correspond to 5, 5' — 100, while 
the dotted curves correspond to S,S' = 0.1. Anti-correlation 
is clearly exhibited in all cases. 



dressed energies in the infinite-time spectral distribution 
of the output beam. The same kind of coupling is present 
in the bosonic case, albeit its effects are less dramatic if 
only the ground state is initially populated (which is the 
situation when the temperature is zero). 

A peculiar characteristic of the system here consid- 
ered, also present in one-dimensional single-mode boson 
models, is the presence of a bound mode of the coupled 
system, for any value of the coupling constant. This im- 
plies a non-Markovian behavior of this system, and has 
two important consequences: the probability of finding 
the atoms in the trap is oscillatory, and a fraction of the 
atoms remains in the trap, even in the infinite-time limit. 
For strong coupling, this fraction approaches 1/4 when 
the number of atoms is much larger than one. This bound 
mode should be however highly sensitive to an external 
potential like a gravitational field. 

Since in the model here considered there is no external 



replenishing of the trap, the outgoing beam has a non- 
stationary nature. For finite times, it displays a very 
intricate dynamics, which results from the combined ef- 
fect of a train of wave packets, with transient behavior 
that exhibits the so-called "diffraction in time" effect, 
and which overlap with the bound-state wave function. 

In the infinite-time limit, however, it is possible get 
an analytical expression for the spectral distribution of 
the outgoing atoms. In the strong-coupling limit, and 
for a steep trapping potential, the outgoing atomic beam 
exhibits remarkable features, for large times: it is quasi- 
monochromatic, and it displays anti-bunching. It is in- 
teresting to remark that the combination of these two 
features is highly desirable, although hard to achieve, in 
light beams. Indeed, the generation of low-noise laser 
light has been an intense field of research [sj , since the 
first experimental observations of anti-bunching [s^ and 
sub-Poissonian statistics 's^. For the fermionic beams 
considered here, anti-bunching comes out quite naturally. 
On the other hand, we have shown that, under certain 
conditions, it is also possible to get here, in the infinite- 
time limit, a quasi- monochromatic spectral distribution, 
in spite of the large number of occupied energy levels in 
the trap. This could be especially helpful for some ap- 
plications recently envisaged for fermionic atomic beams, 
like for instance the development of low-noise atomic in- 
terferometers [5^. 

The investigation of more realistic situations, including 
for instance the presence of a gravitational field, will be 
the object of further consideration. 
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APPENDIX A: SUM FOR THE INFINITE BOX 

In this appendix, we calculate the sum in Eq. (|^ . 
From Ref. [s^ we have the following result: 



1 °° 1 
7rcot(7rz) = — h ^ 1^ 



1 



so that 



z ^ — ' Lz — m z + m 

m— 1 



00 ^ 
TTCOt(TTz)~Z > — 



(Al) 



(A2) 
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In our case, we have the sum: 



^= E ■ 

m=i(odd) 



(A3) 



which may be written as 

oo ^ 

^ jv>-) 2 



E 



1 



(A4) 



m— 1 



m=i(even) 

Setting in the second sum m = 2n, we may write: 



S ■ 



2! E z2 - m2 ~ E ^5—7^ 



2 L — ' — 

m— — CO n— — 00 



z2 - (2n)2 

so that, from Eq. (|A2|) . 

"cot(7rz) cot(7rz/2) 



Since 



cot (a;) = - 



cot 



(f) 



2z 



X 

— tan I — 
2 



we finally obtain: 

CSO 

E 

m=i (odd) 



TT tan(7rz/2) 



(A5) 



(A6) 



(A7) 



(A8) 



which leads to Eq. (|71|l . 

Also, letting z iz, we get: 



E 

=1 (odd) 



1 TT tanh(7rz/2) 

z2 + Tin? 4 z 



(A9) 



which when applied to Eq. H39|l leads to Eq. H4U|I . 



APPENDIX B: SUM FOR THE HARMONIC 
OSCILLATOR 

In this Appendix, we evaluate the sum in Eq. H4()|) . We 
start by proving the identity 



^-E 



^-'o z + m ^ r(z+l/2)' 



where T{z) is the Gamma function, d — ^Jh/mu)Q is the 
width of the ground state of the harmonic oscillator, and, 
from Eq. l|l5|l , 



</'2m(0) = 



1 



,1/4 



T(P' ^22"(2to)! 



H2m(0), (B2) 



where Hjjx) is the Hermite polynomial of order n, with 
(see ref. 0, p. 777) 



H2m{0) - (-1)' 



(B3) 



Replacing these two last expressions into Eq. I)B1|I , we 
obtain 



00 

^E 



(2to)! 1 



/— « 22™(m!)2 z + m 

m— 



(B4) 



Let us consider now the function M(z) = 
V^r(z)/r(z + 1/2), and prove that S = M{z). The 
function M{z) is a meromorphic function, with poles on 
the non-positive integers in the complex plane: z = — m, 
TO = 0, 1, 2, 3, . . .. Therefore, we may write ,56j : 



, [Residue of M(z) for z = —ml . , 

^(^) = E ttt:^ ■ (B5) 



m=0 



z + m 



The residues in the above equation are given by: 

V^r(z) 



Res[M (z), — m] = lim {z + m) 



z — — rn 



r(z + i/2) 



TO!r(-m + 1/2) 

We use now that: 

1 sin(|-TO7r)r(m + l/2) 

r(l/2-TO) ~ n 



(B6) 



r(TO + 1/2) 



(2m)!r(l/2) 
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(B7) 



(B8) 



and r(l/2) = 0?, and replace Eqs. (|B2l, and ||B8)) 

into Eq. I|B5|) . obtaining finally: 



V^r(z) 



r(z + l/2) 22™(m!)2 z + m 



(2to)! 1 



which proves the desired identity. 
It follows then immediately that 
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